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In light of new data on neutron distributions from experiments with antiprotonic atoms [ Trzcinska 
et at, Phys. Rev. Lett. 87, 082501 (2001)], we reexamine the role of nuclear-structure uncertainties 
■ in the interpretation of measurements of parity violation in atoms using chains of isotopes of the 

same element. With these new nuclear data, we find an improvement in the sensitivity of isotopic 
chain measurements to "new physics" beyond the standard model. We compare possible constraints 
on "new physics" with the most accurate to date single-isotope probe of parity violation in the Cs 
atom. We conclude that presently isotopic chain experiments employing atoms with nuclear charges 
£vQ ' Z < 50 may result in more accurate tests of the weak interaction. 
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I. INTRODUCTION 



Atomic parity non-conservation (lj ||] (PNC) provides powerful constraints on extensions to the standard model of 
c/2 ■ elementary particles in the low-energy electroweak sector. For example, a deviation in the observed weak charge of 
an atomic nucleus, Qw, from the prediction of the standard model may hint at the existence of extra neutral-gauge 
Z-boson. Other possible "new physics" scenarios are discussed, e.g., in Ref. ||. 

The most accurate to date measurement of atomic PNC has been carried out by Wieman and co-workers (J, [f| 



using a single isotope of atomic cesium, 133 Cs. In such measurements one determines a parity- violating signal £p 



PNC, 

related to the weak charge as £pnc = ^pnc Qw- The parameter fcpNC is supplied from sophisticated atomic-structure 
calculations. Even for the relatively well- understood univalent Cs atom, the accuracy of the calculation of fcpNC 
remains the limiting factor in the determination of the weak char ge. 

There is an ongoing discussion in the literature (5[ |[ @, ||, |§, [H], 0, Q about whether the 133 Cs weak charge 
deviates from the prediction of the standard model. This possible deviation may be interpreted as an indication for 
an extra neutral-gauge boson Z' |l3|, [l4| ; there were numerous discussions in the literature about implications of the 
£N) ■ possible deviation. It is clear that independent tests of parity violation in atoms are required at least at the level of 
t— I ' the present 1% accuracy for Cs. It is worth mentioning that PNC measurements were also carried out in Tl jL5l |T(i| . 
Pb |L7|, and Bi (see, e. g., 113)) ■ Among these, the simplest atom is Tl, but even for Tl the theoretical uncertainty 
for fcpNC is a factor of a few larger than that for Cs [|(| . 
C/3 | An alternative approach allowing one to circumvent the difficulties of atomic-structure calculations was proposed 
by Dzuba et al. The main idea was to form a ratio 1Z of PNC amplitudes for two isotopes of the same element, 
thus cancelling out associated uncertainties of the atomic theory. However, Fortson et al. (2^] pointed out a conceptual 
limitation of this method - an enhanced sensitivity of possible constraints on "new physics" to uncertainties in the 
neutron distributions. As an example, the differences between neutron and proton root-mean-square radii for 133 Cs 
differ by a factor of four in relativistic and non-relativistic nuclear-structure calculations and depend on nuclear 
models. Unfortunately, at the present level of theoretical understanding of neutron distributions, such large nuclear- 
structure uncertainties would preclude an extraction of useful information on weak interactions from isotopic ratios 
measured for heavy atoms. 

Given the inadequate accuracy of nuclear-structure calculations for the analysis of PNC measurements based on 
isotopic ratios, here we investigate the role of uncertainties in neutron distributions using empirical data. Recently, 
Trzcinska et al. [ p3| deduced differences between root-mean-square radii R n and R p of neutron and proton distributions 
from experiments with antiprotonic atoms. A wide range of stable nuclei were investigated and the differences were 
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approximated by a linear dependence suggested in Ref. ]24| 

N — Z 

AR np = (-0.04 ±0.03) + (1.01 ±0.15)— fm. (1) 

iV ± Z 

Here AR np = i?„ — Z is the nuclear charge, and N is the number of neutrons. Recently, this result was employed 



to estimate the nuclear-structure uncertainty for parity- violating amplitude in Cs |10 . In light of the new nuclear 
data we reexamine the suitability of isotopic chain measurements for studies of parity violation in atoms. We find 
that the nuclear-structure uncertainty in possible probes of "new physics" with isotopic chains is reduced by the new 
antiprotonic-atom data. We compare constraints on the direct "new physics" with what is currently the most accurate 
single-isotope probe of parity violation in 133 Cs. We conclude that presently isotopic chain experiments with atoms 
having Z < 50 may be competitive with this single-isotope determination. 



II. BACKGROUND 



In a typical atomic PNC setup, one considers a transition between two atomic states \i) and |/) of the same 
nominal parity. The weak interaction admixes the states \n) of the opposite parity, leading to the otherwise forbidden 
parity-violating amplitude 



E 



PNC 



E 



(f\D z \n)(n\H w \i) 



Ei — E n 
(f\H w \n){n\D z \i) 



Ef — E n 



(2) 



where D is the electric-dipole operator and -ffyv is the Hamiltonian of the electron-nucleus weak interaction. As 
demonstrated by Pollock et al. [|25| , matrix elements of H-w m ay be represented as 



(j\Hw\i) 



Gf 
2Vl 



Cji R^- 2 Q W (N,Z), 



(3) 



where factor Cji depends on atomic wavefunctions and 7 = y/l — (aZ) 2 . 

Including the dependence on nuclear shapes, the nuclear weak charge Qw(N,Z) may be represented at the tree 
level as 



Qw = ~N q n + Z q p (1 - 4 sm 2 w ) + AQ n 



(4) 



Here sin 2 9w — 0.23117 (16) H and quantities q n and q p , introduced in p2| , depend on neutron and proton distribu- 
tions inside a nucleus. It should be noted that quantities q n and q p are numerically very close to one. For example, 
in the "sharp edge" model of nuclear density distribution [p2| 



'In 



1 - - (aZ) 2 
70 v ' 



(5) 



More sophisticated expressions may be found in Ref. ^6|, but the accuracy of the above formula is sufficient for the 
goals of the present work. We omitted radiative corrections in the definition of the weak charge, Eq. ([|). These 
contributions are important in the studies of "oblique" corrections, discussed, e.g., in Ref. |2^, 27 1. Here, motivated 
by possible deviation of the Cs weak charge from the prediction of standard model, we analyze constraints on direct 
tree-level "new physics" . The term AQ n cw in Eq. (||) characterizes "new physics" . Following Ramsey-Musolf J2f| , we 
represent it as a combination of couplings to up (u) and down (d) quarks, i.e. 



AQ new = (2Z + N)h u + {Z + 2N)h d 



(6) 



Z h p + N h n , 



where h p — 2h u ± hd and h n = 2hd ± h u are couplings to protons and neutrons. Various elementary-particle scenarios 
for these interactions were reviewed in Ref. p8fl. Finally, 



N h + Z h p + N h v 



(7) 
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n = n °< 1+ Nwi + [w h » +h ») w> (10) 



with 

y 

ho = -q n + — q P (1 - Asm 2 w ) » . (8) 

As an outcome of the analysis of PNC experiments one would like to set bounds on "new physics" couplings h p , h n 
or equivalently h u , hd', below we summarize relevant analysis from Refs. 0, ^5], ^8|. The PNC amplitudes Ep^c an d 
E' PNC are measured for two isotopes of the same element with neutron numbers N and N' = N + AN, and the ratio 
is formed 

n = E PNC = Qw_ f Rp\ 27 2 ^ 

EpNC Qw \R'pJ 

Here all quantities with primes are for the isotope with N' neutrons. Using Eq. (Q) one obtains 

Z AN h p ( Z , u \ h' - h 

/ s 27-2 

with 1Zq = ( -of ) 'WW • ^ as * term in the above expression may be safely neglected and we determine a 

contribution of "new physics" 

/*„ ^0 J ZAN' 1 > 

In the absence of new couplings T — 0. It may seem counterintuitive that the isotopic ratios are sensitive to the new 
physics encapsulated in couplings to protons (h p ) instead of those to neutrons (h n ). The dependence on h p may be 
easily demonstrated with an alternative ratio (Qw /N — Q w /N')/(Qw /N + Q' w /N'); this ratio is straightforwardly 
reduced to ZAN/(2NN')h p . 

The constraints on h p , Eq.( [ll]), are affected by (i) the experimental error bar in 1Z and (ii) uncertainties in TZq 
which are induced by insufficient knowledge of nuclear distributions. Explicitly, 

Here Aq n = q n — q' n and we assumed R p « R' p . The radii of proton distributions are known with sufficient accuracy p9[ ] 
and we disregarded associated uncertainties. Finally, 

NW (511 3 ^ 2 SARn \ 



ZAN [Ko 7 Rp 
with AR n — R' n — R n . The above expression is similar to the results of p5|] . 

III. RESULTS AND DISCUSSION 

We assume below that in Eq. ( p^ ) the experimental errors 5TZ may be neglected in comparison to nuclear-structure 
uncertainty. By contrast to the previous discussions |22|, [2£| of atomic parity- violation in isotope chains we employ 
the empirical Eq. p) t o estimate radii of neutron distributions; this relation was deduced from experiments with 
antiprotonic atoms ||23| . To estimate the error bar in the differences AR np — AR' np , we formed all possible isotope 
pairs from the original 21 point data set of Ref. |2j| and obtained with the least-square method AR' np — AR np = (0± 
0.003) + (1.01 ± 0.04) {(N' - Z')/(N' + Z') - (N - Z)/(N + Z)} fm. Instead of a single-parameter fit, we performed 
a two-parameter fit because there is no strong theoretical reason to believe that the difference AR 1 — AR np should 
vanish for two distinct nuclei with the same relative neutron excess, (N' — Z')/(N' + Z') = (N — Z)/(N + Z). Such 
obtained (statistical) error bars are very small. However, given insufficient information on systematic errors in [[231 , 
in our subsequent analysis we retained more conservative uncertainties from Eq. ([!]). Based on error bars in Eq. (|1|) 
we set 



<5Ai?„ w <5 (AR' np - AR np ) 



N' — Z N 

(0.03)'- + { 0.15 ' 



N' + Z N + Z 



1/2 

fm. (14) 
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TABLE I: Contribution of nuclear-structure uncertainty to a constraint on "new physics" SJ- for representative isotope pairs. 



Atom 




Mass numbers A 




8T x 10 3 


Ba (Z = 56) 


130 




138 


6.2 


Sm (Z = 62) 


144 




154 


6.5 


Yb (Z = 70) 


168 




176 


12 


Pb (Z = 82) 


204 




208 


39 


i33 Cs (Z = 55)" 
133 Cs (Z = 55) b 


133 
133 






3.7 
5.9 



a Single-isotope constraint for extra neutral-gauge Z-bosons scenario, Eq. ( |16| 
b Single-isotope constraint for isoscalar scenario, Eq. (J16I). 



The first (isotope-independent) term in this expression dominates for Z > 20; for heavy atoms 5AR n w 0.03 fm. It is 
worth emphasizing that the Eq. (|l|) for differences between neutron and proton r.m.s. radii was obtained in Ref. |23j] 
with data for stable isotopes; it is expected that nonstable isotopes may exhibit anomalous AR np . 

We require the nuclear-structure uncertainty in 8T be lower than the current limits deduced from the most accurate 
to date single-isotope 133 Cs determination. The single-isotope measurements are sensitive to a different combination 
of new h u and hd, u — e and d— e couplings. For illustration we parameterize h u = Xhd- For example, A = arises in 
analyses of extra neutral-gauge Z bosons in Eg theories and A = 1 corresponds to pure isoscalar couplings [ p8| . We 
obtain 



^(133 Cs) = „ . (15) 



6h p ^ SQw N 
~h^~ ~ Qw Z+ $ 

We set S q™ k, 0.01, i.e. the present 1% precision of the determination of the weak charge in 133 Cs, and find 
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8.3 x 10- 
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A = 


Z+2N 
2Z+N 


K 





(16) 

\.22hd 



In our illustrative example, the single-isotope bounds set on "new physics" encapsulated in h p are clearly model- 
dependent. We note that the single-isotope 133 Cs measurement is insensitive to a particular scenario h u w — 1.22/id, 
which may be directly probed by the measurements with chains of isotopes or constrained by other electroweak 
observables. 

Given an experimental precision 51Z/1Z in determination of PNC amplitudes, the uncertainty (^) may be minimized 
by using a pair of isotopes with the maximum possible spread of neutron numbers AN. Based on Eq. ( |l4| ) and (|l3| ) we 
calculated 8 J- for such stable isotope pairs for Ba, Sm, Yb, and Pb. We have chosen these atoms mostly because PNC 
experiments were carried out for them, or at least discussed in the literature |l|, [32| |3^, |34|, [55|. ^From the results 
compiled in Table ||, it is clear that the present nuclear-structure uncertainty still may cloud a competitive extraction 
of "new physics" from isotopic chain experiments for these atoms. Compared to single-isotope 133 Cs determination, 
measurements with isotopes of Ba and Sm would be two times less sensitive to extra neutral-gauge Z-bosons and 
would have a comparable sensitivity to new isoscalar physics (see Eq. (|l6|).) Possible constraints from heavier Yb and 
Pb would be affected by the nuclear uncertainty to a larger extent. 

Now we proceed with a search for atoms suitable for isotopic ratio experiments given present nuclear-structure 
uncertainties. In Fig. |l| we summarize results for atoms with nuclear charges 40 < Z < 82. To minimize the effect 
of experimental error S1Z in 8T, the spread in neutron numbers A./V should be chosen as large as possible; we only 
considered atoms with stable isotopes so that AN > 8 ( AN = 4 for Pb). We approximate N w 1.5 Z, R p oc Z 1 / 3 , and 
the error SAR n w 0.03 fm. Thus the nuclear-structure uncertainty in the determination of "new physics" 8T grows as 
Z 8 / 3 , explaining a general trend in Fig. [|. We compare the uncertainties to constraints set by the Cs determination 
(horizontal lines). We conclude that the isotopic chain measurements in atoms with nuclear charges Z < 50 may 
provide comparable limits on couplings for the interesting extra Z scenario. For these elements an interpretation of 
the measurements in terms of direct new physics may be relatively free of nuclear-structure uncertainties. It is worth 
emphasizing that extra Z' were discussed recently in connection with a possible deviation of 133 Cs weak charge from 
the prediction of the standard model. 

We would like to briefly comment on the required experimental accuracy in determination of ratio 1Z of the parity- 
violating amplitudes. Approximating N w 1.5 Z we find 
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FIG. 1: Contribution of nuclear-structure uncertainty to possible constraints on "new physics" ST for atoms with nuclear 
charges 40 < Z < 82. Horizontal lines represent limits derived from single-isotope 133 Cs PNC analysis in the isoscalar (dashed 
line) and extra neutral-gauge boson Z' (solid line) scenarios. 



We set 8T to constraints derived from the determination of 133 Cs weak charge, Eq. (|l6|). We arrive at 

STZ / AN , AN 

— < 0.4— 8T{ 133 Qs) « 0.02 — . (17) 

The required accuracy in the ratio of PNC amplitudes 1Z is in the order of 0.3% for Ba and Sm, 0.2% for Yb, and 
0.1% for Pb. The required experimental error is less demanding for lighter atoms. 

So far the most accurate measurement of parity-violating amplitude was carried out in Cs [Q ; the achieved accuracy 
was 0.35%. As first noted by Bouchiat and Bouchiat [3^| , the matrix elements of the weak interaction scale as Z 3 ; 
the parity-violating amplitude may be weaker for atoms with nuclear charges Z < 50 which are lighter than Cs 
(Z = 55). However, the required experimental error in ratios of PNC amplitudes, Eq. (fl7|), is less demanding for 
lighter atoms. Also an enhancement of PNC amplitude may arise due to an admixture to the initial/final atomic state 
of an energetically close intermediate state of an opposite parity by the weak interaction. For example, calculations 
p7| , |38[ demonstrated that the PNC amplitude for the 6s 21 So — > 5d6s 3 Di transition in Yb is approximately 100 
times larger than in Cs. 

We conclude that at the present level of understanding of neutron distributions, atoms with nuclear charges Z < 50 
may be suitable for competitive tests of parity violation with isotopic ratios. If parity- violating enhancement scenarios 
would be realized for such atoms, the experiments may become feasible. It is worth carrying out a systematic search 
for enhanced PNC amplitudes for atoms and ions with Z < 50. Such an atomic-structure search is certainl y a 
nontrivial task, requiring in most of the cases an accurate account of correlations. For example, Xiaxing et al. |3£| 
argued that the PNC amplitude for the 6s 2 — > 5d7s 3 Di transition in Ba is an order of magnitude larger than in 
Cs. Their semi-empirical calculation was based on assumption that an intermediate state 6s7p 1 P°, which is only 258 
cm -1 deeper than hd7s 3 D\ state, provides the main contribution to the PNC amplitude. To verify their conclusion, 
we have carried out the accurate calculation of this amplitude with combined method of configuration interaction and 
many-body perturbation theory E(J. Our determination resulted in the PNC amplitude 30 times smaller than the 
prediction |39| . The main reason for the discrepancy is the strong interaction of the configurations forming 5d7s 3 Di 
and 6s7p 1 P° states which was not accounted for in Ref. p9[ . This configuration interaction leads to significant 
cancellations of different contributions to the matrix element (5d7s 3 Di\Hw\6s7p 1 P°) and decreases the contribution 
to the PNC amplitude by an order of magnitude. 

To reiterate, with the new data from experiments with antiprotonic atoms [ p3| we reevaluated the role of nuclear- 
structure uncertainties in the interpretation of atomic parity violation with chains of isotopes of the same element. 
We find that the nuclear-structure uncertainty is reduced by these new data. We compared possible constraints on the 
direct "new physics" with the most accurate to date single-isotope probe of parity violation in Cs atom. We conclude 
that presently isotopic chain experiments with atoms having Z < 50 may be competitive with this single-isotope 
determination. As the neutron distribution measurements become more refined (see, e.g., Ref. [[Il]) we expect that 
competitive probes of parity violation with isotopic ratios of the same element may become feasible for heavier atoms. 
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